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Abstract of DE1 9902326 

Acoustic signals of rotating machine as well as rotational speed of reference component of the 
machine are recorded. Rotational speed of drive shaft is digitized and analyzed, using known order 
analysis to determine spectral components of continuously occurring acoustic signals, and suitable 
algorithms are used to provide true-to-expectation value estimate of order spectrum. Various status 
changes are measured and compared to a threshold. Operating statuses to be distinguished are 
connected at fixed rotation speed ranges. Maximum rotation speed range is divided into intervals of 
pertinent adaptation and monitoring phases, with pertinent thresholds determined for every operating 
status of machine and change measure. 
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(S) Verfahren zur Schadensfruherkennung von rotierenden Maschinen 

(57) Verfahren zur Schadensfruherkennung von rotierenden Au 
Maschinen durch Auswertung ihres emittierten Schalls 
mit folgenden Verfahrensschritten: 

- Aufzeichnung und Digitaiisterung von akustischen Si- 
gnalen der rotierenden Maschine, 

- luckenlose Bestimmung der spektralen Komponenten 
der akustischen Signale durch Verwendung der Ord- 
nungsanalyse, 

- Berechnung einer KenngroRe A aus den zentrierten und 
gewtchteten Ordnungspegeln Gy, 

- Ermittlung eines varianzgewichteten und durch vorbe- 
stimmte Signalverarbeitungsalgorithmen geglatteten 
Ordnungssummenpegels Ag aus den Kenngrofcen A, 

- Errechnung eines AnderungsmaftesX aus dem Verhalt- 
nis von Ordnungssummenpegeln Ag, welche zu aufeinan- 
derfolgenden Zeitpunkten ermittelt wurden, 

- vergleichen des Anderungsmafces X mit einem vorzu- 
gebenden Schwellwert. 



Ausgewahlte OrdnungspegetveriSule Qber der Zeit 



^^^^ 




100 150 200 250 300 350 400 450 500 
Zeit 



w iMnp<?nRi ipicfrfi n** rv* 9hq iQn«iR/Q 



A 



DE 199 02 326 C 2 



Beschreibung 

[0001] Rotierende Maschinen im Sinne dieser Beschreibung sind Maschinen, die neben fest stehenden Teilen auch ro- 
tierende Teile wie z. B. Zahnrader, Lager und Weilen enthalten. Besondere Bedeutung haben hierbei Fahrzeuggetriebe, 
5 Elektro- und Verbrennungsmotoren, Turbinen sowie Priifstande fur diese Aggregate. 

[0002] Rotierende Teile emittieren Schwingungen und Gerausche. Die geeignete Auswertung dieser akustischen MeB- 
groBen liefert Informationen zum Zustand des Priiflings. 

[0003] Die vorliegende Erfindung bezieht sich auf ein Verfahren zur Erkennung von signifikanten, akut auftretenden 

Veranderungen im akustischen Vernal ten des Priifaufbaus bestehend aus Priif stand und Priifling, welche auf ein aufgetre- 
10 tenes Schadensereignis riickschlieBen lassen. 

[0004] Es sind verschiedene Verfahren zur Maschineniiberwachung und Schadenserkennung wie Risse, Abnutzung, 

Materialermudung in einem friihen Stadium bekannt. Hierzu werden maschinenspezifische MeBgroBen wie z. B. 

Schwingungen und Gerausche ausgewertet, sofern sich der entstehende Schaden in diesen MeBgroBen abbildet. 

[0005] Aus DE 38 12 474 Al ist eine Einrichtung zur Uberwachung von Walzlagern bekannt, die durch eine fehlerspe- 
15 zifische Frequenzanalyse der akustische Signale des Walzlagers und dem Vergleich mit einem korrespondierenden Re- 

ferenzwert die Schadensbeurteilung durchfiihrt. 

[0006] Ein Verfahren unter Verzicht auf eine Referenz wird in DE 32 45 505 Al vorgeschiagen. Hier werden Schallsi- 
gnale im hoch- bzw. tieffrequenten Bereich mittels zwei entsprechend angeordneten Sensoren erfaBt und zu einer GroBe 
verarbeitet, deren Anstieg zu uberwachen ist. GemaB DE ist dieses Verfahren zur Uberwachung von Walzlagern, Getrie- 

20 ben und Turbinenaniagen vorgesehen. 

[0007] Aus DE 195 45 008 Al sind ein Verfahren und eine Vorrichtung bekannt, die eine Uberwachung von periodisch 
arbeitenden Maschinen zur Fruherkennung von Maschinen veranderungen durch permanenten Vergleich aktueller stati- 
stischer Kennwerte mit Referenzwerten ermoglichen. Die wesentlichen Unterschiede dieser DE mit dem hier vorge- 
schiagenen Verfahren liegen in den eingesetzten mathematischen Analysemethoden, der Bildung eines gleitenden Refe- 

25 renzwertes (im vorgeschlagenen Verfahren) zur Eliminierung des Einflusses zulassiger Zustandsanderungen und der feh- 
lenden Moglichkeit (in der genannten DE) zur Erhohung der Sensitivitat des Verfahrens durch Ausnutzung der a-priori 
Kenntnissen uber mogliche Betriebszustande. 

[0008] Im Gegensatz zu DE 195 45 008 Al werden im hier beschriebenen Verfahren statistische KenngroBen, welche 
aus dem Ordnungsspektrum berechnet werden, nicht direkt zur Klassifikation herangezogen. Vielrnehr werden die spek- 
30 tralen Signalanteile mit Hilfe von statistischen KenngroBen konditioniert. Diese Konditionierung hat zur Folge, daB das 
Verfahren selektiv gegeniiber sehr kleinen Anderungen des Priifobjekts wird, welche nicht dem normalen Verhalten ent- 
sprechen. Anderungen des Priifobjekts, die fur das normale Verhalten des Priiflings charakteristisch sind, werden hinge- 
gen unterdriickt. Die normalen Anderungen durfen dabei wesentlich groBere Werte annehmen als die zu detektierenden 
Anderungen. 

35 [0009] Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren der eingangs genannten Art zu entwickeln, das uni- 
versal einsetzbar ist und ohne die Bereitstellung von Referenzwerten arbeiten kann. Das Verfahren kann somit zur Uber- 
wachung von Priiflingen herangezogen werden, uber deren normales Verhalten zum Beginn der Uberwachung noch 
keine gesicherten Kenntnisse vorliegen. Weiterhin konnen Fehlerzustande detektiert werden, die vormals noch nie auf- 
getreten waren, sofern diese eine Anderung der verwendeten MeBgroBen bewirken. 

40 [0010] ErfindungsgemaB wird die gestellte Aufgabe durch die im Anspruch 1 angegebenen Merkmale gelost. 

[0011] Vorteilhafte Weiterbildungen und Ausgestaltungen des Verfahrens sind in den Unteranspriichen 2 bis 7 angege- 
ben. Die mit der Erfindung erzielbaren Vorteile ergeben sich aus der folgenden Beschreibung. 

[0012] ErfindungsgemaB werden akustische Signale der rotierenden Maschine sowie die Drehzahl einer Bezugskom- 
ponente der Maschine, z. B. die Drehzahl der Antriebswelle, aufgezeichnet, digitalisiert und nach der im folgenden be- 

45 schriebenen Signalverarbeitungsmethode analysiert. 

[0013] Die Bestimmung der spektralen Komponenten der akustischen Signale erfolgt liickenlos durch Verwendung der 
bekannten Ordnungsanalyse. Ein Verfahren zur Ordnungsanalyse, das auf digitalem Resampling aufsetzt, ist in Groppe, 
Jonuscheit, Strama, Thoma: "Ordnungsanalyse 1 ', MeBtechnik und MeBsignalverarbeitung, MESSCOMP 1996, expert 
Verlag, 1996, beschrieben. Zur Bereitstellung einer erwartungswerttreuen Schatzung des Ordnungsspektrums werden 

50 geeignete Algorithmen verwendet. Aus den Ordnungsspektren werden die Ordnungspegel C(u, nta) als Verlauf der ein- 
zelnen Ordnungen uber der Zeit nta bestimmt. Dabei entspricht U der Analysezeitdauer, tiber die ein Ordnungsspektrum 
berechnet wird. Die Variable n beschreibt den Zeitpunkt der Analyse in diskreten Schritten von ta. Der Index u bezeichnet 
die betrachtete Linie des Ordnungsspektrums. 

[0014] Fig. 1 zeigt funf ausgewahlte Ordnungspegelverlaufe, wie sie sich vor einem realen Schadensereignis darge- 
55 stellt haben. 

[0015] Die statistischen KenngroBen li nearer Mittelwert mtfu, nta), quadratischer Mittelwert mq(u, ntj und Varianz 
V(u, nta) jedes einzelnen Ordnungspegeis C(u, nta) werden entsprechend Formel [I] gieitend ermittelt. 

m!(u, ntj = k!mi(u, (n - l)ta) + k 2 C(u, ntj, k t + k 2 = 1 [la] 

60 

m q (u, ntj = k^Cu, (n - l)tj + k 2 C(u, ntj 2 , k v + k 2 = 1 [lb] 

v(u, nta) = mq(u, ntj - mi(u, nta) 2 Pc] 

65 [0016] Der betreffende Ordnungspegel wird zentriert und mit einem aus der aktuellen Varianz ermittelten Faktor in der 
Form entsprechend Formel [II] gewichtet, daB Ordnungspegelverlaufe mit hoher Varianz wenig und solche mit niedriger 
Varianz stark in das Klassifikationsergebnis eingehen. 
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C V ( M ; ntj = (C(m, ntj - mtfu, (n - l)ta))/(v(u, (n - 1)1*) + k), k = 0 . . 1 [II] 

[0017] Einerseits werden durch diese Wichtung Anderungen der Ordnungspegel, welche fur den normalen Zustand des 
Prufobjekts charakteristisch sind, unterdriickt. Andererseits werden selbst minimale Pegelanderungen, welche fur den 
normalen Lauf untypisch sind, "verstarkt" und konnen das Klassifikationsergebnis significant beeinflussen. 
[0018] Aus den zentrierten und gewichteten Ordnungspegeln C v (u, nta) wird mit Hilfe der Vorschrift nach Gleichung 
[HI] eine KenngroBen A(nta) berechnet, welche das Verhalten des Prufobjekts liber der Zeit reprasentieren. 

Mnt a )= f^cMnO ,M,,M 2 =0..O max (III) 
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[0019] Mi, M 2 reprasentieren den Ordnungsspektralbereich, der fiir die Analyse herangezogen wird. 
[0020] Fig. 2 gibt einen varianzgewichteten Ordnungssummenpegel A(nta) wieder, wie er sich vor einem realen Scha- 
densereignis dargestellt hat. 15 
[0021] Diese KenngroBen ermoglichen die Detektion verschiedener Fehlerzustande. Plotzliche Anderungen der Kenn- 
groBen lassen auf akute Anderungen des Priiflings schlieBen. Unter Verwendung von Gleichung [IV] wird aus der Kenn- 
groBe A(nta) ein Merkmal ^(ntg) berechnet, welches die sprunghaften Anderungen des Prufobjekts anzeigt. 

MntJ = A s (nta)/A s ((n - l)tj [W] 20 

[0022] Zur Reduzierung des Einflusses aus Alterung und natiirlichem VerschleiB des Prufobjekts wird der varianzge- 
wichtete Ordnungssummenpegel vor Anwendung der Vorschrift nach Gleichung [TV] mittels digitalem FIR-Filter gemafi 
Gleichung [V] geglattet. 

4X w O=X^(/M("-^f + ,F(/):Filterkoeffizienten, N F :Filtergrad IV) 

[0023] Fig. 3 zeigt das somit ermittelte AnderungsmaB X fur den varianzgewichteten Ordnungssummenpegel, wie er 30 
sich vor einem realen Schadensereignis dargestellt hat. 

[0024] Das vorgeschlagene Verfahren kann fur beliebige rotierende Maschinen ohne jegliche konstruktiven Vorkennt- 
nisse eingesetzt werden. Wird die Uberwachung einer Maschine eingerichtet, so sind Akustik-Sensoren, z. B. Korper- 
schallaufhehmer und Mikrophone, an den Maschinen geeignet anzubringen. 

[0025] Zu Beginn einer Uberwachungsaufgabe werden die akustischen Signale in oben beschriebener Weise analy- 35 
siert. Das AdaptionsmaB 0, welches gemaB [Via] und [VIb] berechnet wird, stellt sich dabei automatisch in einer "Ad- 
aptionsphase" auf einen Bezugswert ein. 

a(nt a )= A(nt a )-A((n-l)t a ) CVIa) 

40 

®0O= +/)«) ,^(/):Filterkoeffizienten, N F rFiltergrad (VIb) 

[0026] Die Dauer der Adaptionsphase ist grundsatzlich abhangig vom zu uberwachenden Maschinentyp, da das An- 45 
laufverhalten der Maschine, bis stationare Betriebsbedingungen erreicht sind, groBen EinfluB auf die Varianz der Ord- 
nungen haben kann. Das vorgeschlagene Verfahren uberwacht die Entwicklung des AdaptionsrnaBes 0, und sobald 0 
iiber eine Zeitdauer T = lta unter einer festzulegenden Schwelle liegt, gilt die Adaptionsphase als beendet und die Uber- 
wachungsphase beginnt automatisch. 

[0027] Uberschreitet das AnderungsmaB X in der Uberwachungsphase einen extern vorzugebenden, auf Erfahrung be- 50 
ruhenden Sen well wert, so ist ein Schadensfall erkannt worden. 

[0028] Uberschreitet das AdaptionsmaB 0 in der Uberwachungsphase einen vorzugebenden, auf Erfahrung beruhen- 
den Schwellwert, so wird X als nicht aussagekraftig bewertet, da das System sich auf das veranderte Verhalten des Priif- 
lings adaptieren muB. 

[0029] Konstruktive Kenntnisse iiber die Maschine konnen genutzt werden, um z. B. Signalbandbreite und Ordnungs- 55 
auflosung fiir eine vorteilhafte meBtechnische Umsetzung des Verfahrens festzusetzen. 

[0030] Eine mehrkanalige Realisierung des Verfahrens kann durch die konstruktive Beschaffenheit des Prufobjekts 
sinnvoll sein, um auftretende Fehler im Sensorsignal besser zu erfassen. 

[0031] Eine vorteilhafte Ausgestaltung des Verfahrens beriicksichtigt die endliche Zahl von Betriebszustanden, in de- 
nen die Uberwachung erfolgt. Das Prinzip dabei ist, daB fur jeden Betriebszustand Bi der Maschine ein AnderungsmaB X{ 60 
bestimmt wird, das beim Erstanlauf des Betriebszustandes Bi in einer Adaptionsphase eingestellt wird. 
[0032] Eine spezielle Ausfuhrung ist, daB die zu unterscheidenden Betriebszustande an feste Drehzahlbereiche gekop- 
pelt sind. Der maximale Drehzahlbereich dmax der zu uberwachenden Maschine wird dabei in I festzulegende Drehzahl- 
intervalle aufgeteilt, wobei die Drehzahlintervalle sich auch uberlappen oder auch luckenhaft aneinandergereiht sein dur- 
fen. Das vorgeschlagene Verfahren ermittelt fur jedes Drehzahlintervall D[ ein AnderungsmaB Xi, das beim Erstanlauf 65 
des Drehzahhntervalls in einer Adaptionsphase eingestellt wird. 

[0033] Der Vorteil dieses Vorgehens liegt in der erhohten Sensitivitat des Verfahrens, da bei der Bestimmung des An- 
derungsmaBes X verschiedene Betriebszustande (und damit unterstellterweise unterschiedliches akustisches Verhalten) 
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unterschieden werden, die nicht mit einer Schadigung der Maschine im Bezug stehen; X[ wird insgesamt "giatteren" Ver- 
lauf haben. Die Adaptionsphasen fur wiederholt auftretende Betriebszustande bzw. Drehzahlintervalle konnen verkiirzt 
werden, wenn das letztmalige Auftreten des betreffenden Betriebszustandes Bi bzw. Drehzahlintervalls Di nicht langer 
als die Zeitdauer Tfii bzw. Toi zuriickliegt und der "alte" varianzgewichtete Ordnungssummenpegel Asi((n - l)ta) noch als 
5 Referenz zur Berechnung von Xi angenommen werden kann. Die Zeitdauer T B i fiir die I Betriebszustande sind durch Er- 
fahrungswerte festzusetzen. Fig, 4 verdeutlicht das Prinzip zur EinteiLung des maximalen Drehzahlbereichs in Drehzahl- 
intervalle und den zugehorigen Adaptions- und Uberwachungsphasen, 

[0034] Eine vorteilhafte Ausgestaltung des Verfahrens beriicksichtigt, daB von bestimmten zu uberwachenden Kompo- 
nenten des Prufobjekts determinierte, a-priori bekannte Ordnungsspektralanteile emittiert werden. Zur sensitiven Uber- 
10 wachung der einzelnen Komponenten einer rotierenden Maschine unabhangig voneinander konnen fiir mehrere Satze 
von Ordnungslinien Mj jeweils ein AnderungsmaB Xj ermittelt werden. Unter Beriicksichtigung verschiedener Betriebs- 
zustande Bi bzw. Drehzahlbereiche Di bedeutet dies, daB in alien Betriebszustanden Bi bzw. Drehzahlbereichen Di fur je- 
des Mj ein AnderungsmaB Ay berechnet wird, das beim Erstanlauf des Betriebszustandes Bi bzw. Drehzahlintervalls Di in 
einer Adaptionsphase eingestellt wird. 

15 

Patentanspriiche 

1. Verfahren zur Schadensfruherkennung von rotierenden Maschinen durch Auswertung ihres emittierten Schalls 
mit folgenden Verfahrensschritten: 

20 - Aufzeichnung und Digitalisierung von akustischen Signalen der rotierenden Maschine, 

- luckenlose Bestimmung der spektralen Komponenten der akustischen Signale durch Verwendung der Ord- 
nungsanalyse, 

- Berechnung einer KenngroBe A aus den zentrierten und gewichteten Ordnungspegeln Cy, 

- Ermittlung eines varianzgewichteten und durch vorbestimmte Signalverarbeitungsalgorithmen geglatteten 
25 Ordnungssummenpegels A« aus den KenngroBen A, 

- Errechnung eines AnderungsmaBes X aus dem Verhaltnis von Ordnungssurnmenpegeln Ag, welche zu auf- 
einanderfolgenden Zeitpunkten ermittelt wurden, 

- vergleichen des AnderungsmaBes X mit einem vorzugebenden Schwellwert. 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB zur sensitiven Uberwachung einzelner Maschinen- 
30 komponenten unabhangig voneinander fur mehrere Satze von Ordnungslinien Mj jeweils ein AnderungsmaB Xj er- 
mittelt wird. 

3. Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, daB durch eine mehrkanalige Anordnung die kon- 
struktiven Eigenschaften der zu uberwachenden Maschine insofern besser beriicksichtigt werden, als daB die Sen- 
soren naher am moglichen Fehlerort sind und damit der akustische EinfluB des Fehlers im Schallsignal verstarkt ent- 

35 halten ist. 

4. Verfahren nach einem der Anspriiche 1 bis 3, dadurch gekennzeichnet, daB fur jeden Betriebszustand Bi, der Ma- 
schine ein AnderungsmaB Xij bestimmt und ein zugehoriger Schwellwert vorgegeben wird. 

5. Verfahren nach Anspruch 4, dadurch gekennzeichnet, daB die Adaptionsphase fur wiederholt erreichte Betriebs- 
zustande Bi entfallt, indem der betriebszustandspezifische, varianzgewichtete Ordnungssummenpegel A^^n - l)ta) 

40 als Referenz zur Berechnung von ^j weiterverwendet wird, sofern das wiederholte Anlaufen des Betriebszustands 

Bi nicht spater als eine vorzugebende Zeitdauer T B i erfolgt. 

6. Verfahren nach einem der Anspriiche 1 bis 3, dadurch gekennzeichnet, daB der maximale Drehzahlbereich der 
Maschine in Drehzahlintervalle unterteilt wird, und fur jedes Drehzahlintervall Di ein AnderungsmaB A*,- bestimmt 
und ein zugehoriger Referenzwert vorgegeben wird. 

45 7. Verfahren nach Anspruch 6, dadurch gekennzeichnet, daB die Adaptionsphase fur wiederholt erreichte Drehzahl- 

intervalle Dj entfallt, indem der betriebszustandspezifische, varianzgewichtete Ordnungssummenpegel Agi((n - l)ta) 
als Referenz zur Berechnung von ^j weiterverwendet wird, sofern das wiederholte Anlaufen des Drehzahlintervalls 
Di nicht spater als eine vorzugebende Zeitdauer T Di erfolgt. 
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Figur 2 
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